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Abstract 

We derive recurrence relations for the calculation of multiloop sunset-type 
diagrams with large powers of massive propagators. The technique is for- 
mulated in configuration space and exploits the explicit form of the massive 
propagator raised to a given power. We write down and evaluate a convenient 
set of basis integrals. The method is well suited for a numerical evaluation 
of this class of diagrams. We give explicit analytical formulae for the basis 
integrals in the asymptotic regime. 
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1 Introduction 



With the present accuracy of phenomenological applications, high precision tests of 
the Standard Model and search for new physics invariably require accounting for 
higher order contributions in coupling constants within perturbation theory This 
necessitates the computation of multiloop diagrams M. Numerous high precision 
tests of the Standard Model require the evaluation of diagrams in the approximation 
when external momenta are small in comparison with the masses of particles corre- 
sponding to internal lines of the diagrams (see e.g. |2], [5], ||). In this way massive 
diagrams without external momenta - vacuum bubbles - appear. 

The dominant technique for high-order perturbation theory calculations of mul- 
tiloop diagrams nowadays consists in using recurrence relations obtained from the 
method of integration by part within dimensional regularization ||. The final an- 
swer of a calculation is expressed in terms of a few master integrals that serve as 
initial values for the recurrence equations. These master integrals are explicitly eval- 
uated (analytically or numerically) while the reduction of a given diagram to the 
master integrals involves only algebraic manipulations (as an example see ref. |J 
which presents the idea in its full bright glory). 

Presently the evaluation of massive vacuum bubbles diagrams at the three-loop 
level is basically completed. The general strategy of reducing all massive three-loop 
bubbles to a basic set of master integrals through recurrence relations was described 
in ref. 0. Recently the remaining unknown master integrals have been analytically 
identified by using high precision numerical computation assuming a given basis 
of transcendental numbers for master diagrams ||. The existing techniques for 
the three-loop evaluation of massive vacuum bubbles make use of the reduction 
of six-dimensional objects (in the parameter space of powers of the propagators) 
through recurrence relations obtained within the integration by part technique. In 
direct applications of these techniques some topology classes of diagrams (sunset- 
type topology with large powers of propagators) appear to be one of the main sources 
of heavy computer time consumption ||. Nevertheless, in principle, the problem of 



2 



analytical calculation of three-loop bubbles is solved. 

However, in practice the solution to the recurrence relations is not a routine 
procedure and requires the manipulation of a huge number of terms. Even using 
symbolic manipulation programs this task is in some cases beyond present computer 
capabilities. For instance, the computation of higher moments of the 6-quark spec- 
tral function fTO , |TT| , which is important for a precision determination of the 6-quark 



mass 



12| , is limited by these obstacles and only moments up to n = 8 are presently 
available within perturbation theory [[HJ. In view of the importance of this problem 
new ideas and techniques to improve on present results are called for. 

Some attention has recently been drawn to the problem of optimizing the recur- 
rence procedure for three-loop bubbles to find shorter routes to the final solution 
(see e.g. |14| and refs. therein). It has been proven that the final result for three- 
loop bubbles contains only several transcendental numbers (which are known) with 
rational coefficients for arbitrary mass configuration ||. This structure of the fi- 
nal answer prompts one to search for more direct ways of obtaining the physical 
results. Indeed, for phenomenological applications one only needs numerical values 
of the coefficients of perturbation theory and their analytical expressions are not 
really important ||. The advantage of analytical computation is, of course, the 
full control over precision. If numerical methods are used, one has to take extra 
care about possible error accumulation. This is not always straightforward in this 
type of calculation because huge numerical cancellations may occur among the large 
number of contributing terms. 

We discuss a technique that allows one to reduce these diagrams to a one- 
parameter set of basis integrals of a rather simple structure. The method is formu- 
lated in configuration space and exploits the explicit form of massive propagators 
with large powers of denominators. A convenient set of basis integrals is written 
down. The generalization to any number of loops is straightforward which is a part 
of the motivation for the present investigation. The method is well suited for a 
numerical evaluation of this class of diagrams. 

In the course of comparing our analytical results with those given in ref. [0] 



we have derived a number of interesting new analytical results on definite integrals 
involving products of McDonald functions with powers and logarithms. 

The paper is organized as follows. To begin with, in Sec. 2 we remind the reader of 
some features of the configuration space approach and introduce our notation which 
closely follows the notation in ref. [JT5| . In Sec. 3 we present several explicit examples 
taken from the subclass of sunset-type diagrams, the water melon diagrams, in order 
to compare the results obtained within the configuration space approach with results 
given in the literature ||. In Sec. 4 we explain our general recursion concept and 
present a convenient set of basis integrals. In Sec. 5 we present different methods to 
calculate these basis integrals. Our conclusions are given in Sec. 6. 



2 Basic relations 

We consider three-loop vacuum bubbles with only one mass m which serves as a 
dimensional parameter. The classification of the topology prototypes for three-loop 
vacuum bubbles was presented in ref. 0. The analytical computation of some 
missing master integrals has recently been completed ||. However, the solution 
of the recurrence relations leading to the master integrals is complicated and time 
consuming, especially for large powers of propagators. We suggest new recurrence 
relations for a particular topology of vacuum bubbles which allows for an explicit 
solution. The simplicity of our technique is manifest in the configuration space 
representation for Feynman diagrams. First we remind the reader of some features 
of the configuration space approach and introduce our notation which closely follows 
the notation in ref. |15] . 



A general three-loop vacuum diagram has the topology of a tetrahedron. It 
also can be thought of as "fish+propagator" topology (see Fig. 1) where the "fish" 
part is a master two-loop diagram. This topology suggests the use of the following 
representation in configuration space (with an obvious choice of space-time points) 

IL(x) = ftsh(x)D(x,m) (1) 
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(a) (b) 

Figure 1: Three-loop vacuum bubble diagram in two different representations: (a) 
the tetrahedron representation, (b) the "fish+propagator" representation where the 
configuration space points and x are indicated. 

where D(x, m) is a propagator of a massive particle with mass m in D-dimensional 
(Euclidean) space-time, 

1 r e ~^d D p _ (mx) x K x (mx) T(\) 
[X,m) (2n) D J p 2 + m 2 {2ir)^x™ ' ] Arr^x^ [ ) 

We write D = 2A + 2, A = f — e; K\(z) is the McDonald function (a modified Bessel 
function of the third kind, see e.g. |L6|, [T?|). The propagator D(x,m) depends only 
on the length of the space-time vector \x\ = y/x^x 11 for which we simply write x. 
The explicit representation for the modified internal line with mass m (power of the 
massive propagator) is given by 

(X ' ^ = J (p2 +m 2)M+l = (27T) A +W(/X+1) \ X~) K ^ mX >- 

(3) 

It contains the same functions (up to the difference in indices) as Eq. (fj) and thus 
does not change the general functional structure of the representation constructed 
below. This is the reason why our method is well suited for dealing with large powers 
of propagators of massive particles. 

The quantity of interest - a vacuum bubble - is represented by the integral 

f[(0) = f Tl(x)d D x (4) 



which is nothing but the Fourier transform of IT(x) at zero momentum. It depends 
on a single dimensional parameter m, the mass of the particles in the massive lines. 

The "fish" part of the diagram is written in terms of a dispersion relation in 
configuration space, 

fish(x) = J pf(s)D(x,\ / s)ds (5) 
which leads to a representation of the form 

f[(0) = J Tl(x)d D x = J dspf(s) J D(x,^)D(x,m)d D x. (6) 

This form can be further simplified by performing the integration in x explicitly. 
The result reads 

'mv / s\ A 1 r°° 



/( 7Tb\J S \ 1 /*oo 

D(x, y/s)D(x,m)d D x = f — — J A+im T fj J xK x {x^)K x {mx)dx 



(2vr) A + 1 r(A + 1) Jo 

= ("I") (2,)^.r(A + 1) \2 ) \2 ~ ) (2 + A ' ? l! 1 - Tj ' 

Here 2-^1(0, b; c; z) is a hypergeometric function. Therefore the problem of evaluating 
the diagram is reduced to a one-dimensional integral if the spectral density for the 
"fish" part Pf(s) is known. The spectral density Pf(s) was computed for the principal 
massive configurations in four- dimensional space-time [IS]. If the propagator that 
multiplies the "fish" part of the diagram is massless, this formula can be further 
simplified and is given by 

The final integration over s in Eq. @ now includes only some power of the energy 
square s (instead of the hypergeometric function as in the massive case) and is 
rather straightforward. The equivalent set of formulae can also be obtained in the 
momentum space representation ||. 

In some cases one propagator can be removed from a diagram using the recur- 
rence relations for bubbles obtained within the integration by part technique 0. 
Then the diagram becomes simpler. The initial "fish+propagator" topology con- 
verts to a "spectacles+propagator" topology. The spectral density for the spectacle 




(a) (b) 

Figure 2: The "spectacle+propagator" representation, also called "spectacle" topol- 
ogy diagram in two different forms, namely (a) the form used in Eq. fl69|) and (b) 
the form used in Eq. (jSll). 

part can be computed in a rather simple way. It is given by a product of two one-loop 
integrals in the momentum space representation (see Fig. 2). In the original classifi- 
cation of these are class E diagrams. We will henceforth adopt the classification 
of to denote the different topology classes of diagrams. 

After a deliberate use of the recurrence relations for bubbles, in some cases 
two propagators can be removed. A typical situation of such a kind was analyzed 
in ref. 0. In such a case the diagrams become simple indeed. Even the most 
complicated ones belong to the subclass of water melon topologies which can be 
computed immediately (see Fig. 3). Their properties (with generalization to any 
number of loops) have been described in detail in ref. ||15| . The most attractive 
feature of such a strategy is that for high derivatives of propagators (large powers 
of denominators) the corresponding recurrence relations for this particular topology 
can be solved very efficiently 

3 Water melon topology: explicit examples 

While the water melon class of diagrams can appear as part of the remnants of the 
general recursive procedure, there are some cases when they are just the final aim 
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(a) (b) 

Figure 3: (a) Three-loop (i.e. four-line) water melon and (b) two-loop water melon 
(three-line water melon, being the ordinary sunset diagram). 



of the recursion procedure. This is the case for the Bn subclass of diagrams 0. 
Some of the master integrals (for instance, D^O, 1, 0, 1, 1, 1) in ref. Jj]]) are exactly 
water melons. A further simplification of water melon diagrams can be achieved 
with the use of their particular properties. In the configuration space representation 
the water melon diagrams can be reduced to a specific basis set of simple integrals 
quite efficiently 

The configuration space technique for the water melon topology is particularly 
convenient for a numerical evaluation since it is simple and efficient. We have 
reproduced numerically some results of refs. || with our technique. 

In the particular case of the subclass Bn of the bubble diagrams, the water melon 
topology diagrams emerge naturally and can be chosen as master configurations. 
The reduction of a general diagram of this subclass to the water melon topology 
is explicitly constructed in ref. ||. Here we discuss a representation which can be 
obtained within the configuration space technique. The representation is simple 
and is given by a one-dimensional integral. The expansion in e within dimensional 
regularization is straightforward and is explicitly given for the evaluation of the 
numerical value of the integral 5 4 || for which we give a new representation. 

The starting point of our calculation is the definition of the Bn class of dia- 



S 



, d D kd D ld D p 

J B A r(0,0,rz3,n 4 ,n 5 ,n 6 ) = / x 



m 3D( 7r D/ 2r ( 1 + £ )) E 

(7) 



((p + A;) 2 + m 2 )™3 ((p + Z) 2 + m 2 )"4 ((p + fc + /)2 + m 2)n 5 ^2 + m 2)n 6 

with two propagators absent (ni = n2 = 0) to obtain a water melon topology 
for the three-loop case. In the following we will suppress these first two indices 
in the notation for the class diagrams. The configuration space expression 
for the generalized propagator (with crosses or having differentiated in its mass or 
momentum) is given by Eq. (|3]). This can be inserted into the above expression for 
-Btv and after rearrangement of integrations leads to 

™2(n3+n4+n 5 +rt6)— 3D 

B N (n 3 ,n^,n 5 ,n 6 ) = + g))3 X (8) 

(2tt) 3D J £>("3-i) ( X) m )£)(«*-i) ( X) m)^)^- 1 ) ( x , m)/^™ 6 " 1 ) (x, m)d D x 

which can be reduced to a one- dimensional integral using the rotational invariance 
of the integration measure in Euclidean space-time, 

d D x = -x 2X+1 dx. (9) 

T(A + 1) V ; 

Note that various techniques of eliminating tensorial structures for vacuum diagrams 



were discussed in refs. |T5|, |19]] . On the other hand we have [jig] 



2ir x+1 r°° 

11(0) = — / D ( - n3 - 1 \x,m)D^- 1 \x,m) x 



T(A + 1) Jo 

D^- 1 \x,m)D^- 1 \x,m)x 2X+1 dx. (10) 

The comparison of these two formulas results in 

B N (n 3 ,n 4 ,n 5 ,n 6 ) = ^d/2 T{1 + e))3 S (0) (11) 

where the powers of the propagators in 11(0) have been appropriately adjusted. In 
the following calculations we set m = 1 for the mass (regardless of the units, of 
course). We will discuss some explicit applications of Eqs. flip] ) and fllTD in the 
following to indicate their properties. 



As a first example we consider the integral -Bat (2, 2, 2, 2) in the case A = 1 (four- 
dimensional space-time). Generally we have the explicit expressions 



x 



l-A 



(2 7 r) A + 1 2 1 r(2) 
for the propagators which results in 

D^(xA) = 



K x .Ax) 



2(2vr) 



;K {x) 



for e = 0. We obtain 



5^(2,2,2,2) 



(2tt 



,12 



2tt 2 



7T 



16(2tt) 8 Jo 



Kn(x)x 3 dx = 2 Kn(x)x 3 dx 



(12) 



(13) 



(14) 



which is the expression obtained earlier in f!5 |. Note that the function K (x) is 
a propagator of a massive particle in two-dimensional space-time. Therefore many 
results can be obtained by using two-dimensional field theory in Euclidean space- 



time (see e.g. ref. [^0 



For the more general case of D- dimensional space-time we obtain 

2tt 2 - £ 



n(0) = , — / x 4e K\(x)x 3 - 2£ dx 

v ; 16(27r) 8 - 4e r(2 - e) Jo ~ eK ' 



and 



5^(2,2,2,2) 



)l-2e 



K_ £ (x)x 



3+2e 



dx. 



(15) 



'l-s)Y{l + efT{l-e) Jo 
To find higher orders in the e-expansion necessary for computations within dimen- 
sional regularization we use series expansions in e of all quantities entering Eq. (|i~5|) 

First we have the rather obvious results 
2 i- 2e 



;i- e )r(i + e) 8 r(i- e ) 

x 3+2e 



2(1 + e - 2e In 2 + 2e-f E ) + 0(e' 
x 3 (l + 2e\nx) + 0(e 2 ), 



(16) 



where 7^ is Euler's constant. Within the dimensional regularization scheme the 
propagator in the configuration space contains the McDonald function with a non- 
integer index depending on the regularization parameter e. To expand the McDonald 



function in the parameter e entering its index, we use the general formula |T7| 



dK v {z) 
dv 



ne{0,l,...} (17) 
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for the derivative of the McDonald function with respect to its index near integer 
values of this index. In this case one obtains 



K„ £ (x) = K (x) + O(e 2 ). 

We end up with 



yoo poo 

5^(2,2,2,2) = 2/ K*{x)x 3 dx + 2e{l + 2\n2 + 2 lE ) / K*(x)x 3 dx 
Jo Jo 

+4e / K^x)x 3 \nxdx + 0{s 2 ) 
Jo 

= 2I (3) + 2e(l-2\n2 + 2 lE )I (3)+4eI l (3) + O(e 2 ) (19) 

where a general notation for the configuration space integrals 

I m (q) = / Kt m (x)Kr(x)x^dx, 
Jo 

I l m (q) = / K*- m (x)K™(x)x q \nx dx (20) 
Jo 

is introduced. The only new contribution in the ^-expansion till the first order in 
Eq. (|i~9"D is connected with the logarithmic integral I l (3)(P m (q) from Eq. (|20|)). The 
part related to /o(3) in this order is a trivial kinematic contribution. The term 
2(ln2 — 7#) in Eq. (plj) can be easily removed by redefining the logarithmic integral 
using lnrr — > ln(xe 7JS /2). 

Identifying the parameters B 3 and 5 4 from || we find 

5^(2,2,2,2) = -^ + ^3+(|i?3 + ^4)5 + 0(e 2 ) (21) 

= _3 + ! c(3) + ^ c(3) _ 63 c(4) + ± B \ £ + 0(e 2) 
8 16 SV ' V32 sv ' 32 sv 7 16 7 V ' 

where B 3 = ((3) — §sC(4) + 0(e 2 ). The comparison of the zeroth order term of 
Eq. ( |2~T| ) with Eq. (|TJ|) results in the relation 

'o(3) = ~ + !«3) (22) 

which assigns a value to one of the initial terms of the recurrence relations that will 
be presented later in Sec. 4. We checked this relation numerically. In the first order 
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part of the ^-expansion we solve for B4 obtaining the representation 

B, = y(2(l-21n2 + 2 7s )/ (3)+44(3) + g(C(4)-C(3)) 

= |((l-21n2 + 2 7E )/ (3) + 24(3)) + ^(C(4)-C(3)) (23) 

which after substituting for I l (3) from Eq. (p0|) gives numerically _B 4 = —1.7628 . . . 
This numerical value expressed in terms of configuration space integrals within our 
technique coincides with the result given in |j) . Taking the analytical expression for 
B 4 from || 

B, = 16Li 4 Q + I ln 4 (2) - ^ 2 ln 2 (2) - ^ (24) 
with Li 4 (z) being a fourth order polylogarithm, 

00 z k 

k=i K 

we obtain the result for the logarithmic integral /q(3), 

~ poo 

4(3) = / K*(x)\n(xe~< E /2)x 3 dx (25) 
Jo 

= A + -Li 4 (-) - - -7r 2 ln 2 (2) + - ln 4 (2) + —((3) 

32 4 \2J 1920 32 w 32 w 128 v ; 

which serves as the initial value for the recurrence relations for the set of logarithmic 
integrals (Sec. 4). We have checked this value numerically. For a further demon- 
stration of the efficiency of the configuration space technique for fixed powers of 
the propagators we calculate the two integrals 5^(2,2,2, 1) and B N (2, 3, 3, 4) (the 
latter does not contain I? 4 which is the reason for having selected this example). 

For the integral B^(2, 2, 2, 1) in four- dimensional space-time (A = 1) we have to 
include the propagator 

D<0> <*< = (aryCrp) *^ = (f)^ A '- W (26) 

equal to Ki(x)/47r 2 x for A = 1 which is a standard propagator of a massive particle 
for D = 4. We obtain a representation of the form 

B N (2, 2, 2, 1) =4 J Kl(x)K l {x)x 2 dx = 4/ 4 (2). (27) 
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For the integral B^(2, 2, 2, 1) in the case A = 1 — e we obtain the generalization of 
Eq. (|27|) of the form 

B ~< 2 - 2 - 2 - 1] = (i- £) r(i + E )»r(i- E ) 1 K *-W jr '- ix)x2+2 ' dx - (28) 

The e-expansion of the factor multiplying the integral is the same as in the former 
case except for an overall factor of 2. A similar statement is valid for the expansion 
of the power of x. What remains to be done is to expand the McDonald functions 
in the vicinity of integer values of their indices. To obtain this expansion we use the 
relation 



dKJx) 



i rx\^ „ , s i 



dv 

which contributes to the power expansion as 



2V2 / Ko(x) = -K Q (x) (29) 



K x ^{x) = K x {x) - -K {x) + 0{e 2 ). (30) 



x 



Using these expansions we obtain the representation in terms of our basic integrals 

poo 

5^(2,2,2,1) = 4/ x 2 Kl{x)K 1 (x)dx 
Jo 

POD 

+4e(l-21n2 + 2 7£ ) / x 2 K$(x)K 1 (x)dx 

Jo 

POO poo 

-As xK^(x)dx + 8e x 2 \nx K^(x)K 1 (x)dx + 0(e 2 ) 
Jo Jo 

= Ah{2) + 4e(l - 2 In 2 + 2 7b )/ 1 (2) - 4eJ (l) + &el[(2) + 0(e 2 ) 

that has to be compared with the output of the RECURSOR package || with the 
explicit master integrals B 3 , 5 4 , 

5^(2,2,2,1) = 7 -B 3 + ^B 4 e + 0(e 2 ) 

= Ic(3)+^i?4-?C(4)W^ 2 ). (31) 



4 av ' V4 

The zeroth order comparison gives the result ii(2) = ^C(3) which has been verified 
numerically. The first order comparison results in 

A f\ Q 

5 4 = -(4(l-21n2 + 2 7£ )/ 1 (2)-4/ (l) + 8/|(2) + -C(4)) (32) 

16 21 
= _((i_21n2 + 2 7js )J 1 (2)- J (l) + 27{(2)) + -C(4) = -1.7628 
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as before (see Eq. ( p3|) ). For the integral B N (2, 3, 3, 4) in the case A = 1 — e we need 
to include further propagators. They are 

D {2) (x,l) = ——^-—K x _ 2 (x) = -^——K^Ax), 
y ' > (27r) A+1 2 2 r(3) V ; 8(2tt) 2 - £ eK h 

z^>3 A | s 

v ' ; (27r) A+1 2 3 r(4) 3V ; 48(2vr) 2 - £ eK ' 
Both McDonald functions have to be expanded in their index. We have K_\{x) = 
Ki(x) and 

dKAx) 1 fxy 1 j, , s 1„ , s 



thus 



and 



K^ £ (x) = K^x) + -K (x) + 0(s 2 

x 



u=-2 2 1 2! \2J v ; 1!1! \2J v ' 



dv 



2 2 

— riCo(^) Ki(x) 

X A X 



so that 

—K x {x) + -= 
X x z 

With these relations we obtain the following ^-expansion for the integral in question 



K- 2 - 6 = K 2 (x) + -K^x) + -K {x) + 0{e 2 ). (33) 



2-2e 

7^(2,3,3,4) = ; r— r x 

V ' ' ' ; 192(1 -e)r(l + £) 3 r(l -e) 



POO 

/ ir_ £ (x)K! 1 _ e (x)ir_ 2 _ £ (x)a; 7+2£ rfa; 

JO 

1 :hi(7) + -|-(1-21ii2 + 2 7b )/ 21 (7) + -§7/11(6) 



192 v ' 192 v ' ' y ' 96 

+ s /3(6) + s /2(5) + s 4(7)+ ° (e2) (34) 



where we have introduced the generalized integral 



IX, 



I mn (q) = / r""^)^^)^ {x)x«dx. (35) 

This expansion has to be compared with the representation through master integrals 
found in momentum space, 

B „ (2 , 3 ,3,4)^ + (J|_ B3 -^) £ + ( a (36) 
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resulting in the identification 

poo I 

/ 2 i(7) = / K Q {x)Kl{x)K 2 {x)x 7 dx = - (37) 
Jo 3 

which is surprisingly simple and contains no transcendental numbers usually present 
in such integrals. It is rather curious that a similar identification allows one to 
express £(3) in terms of our basis integrals 

C(3) = ^((l-21n2 + 2 7s )/ 2 i(7) + 2/ 11 (6) 

+2/3(6) + 2/ 2 (5) + 24(7) + — ). (38) 

We checked both equations numerically to make certain that they are valid. These 
results serve as a hint that the standard basis may not be the simplest and most 
relevant basis to be used in computations of massive three-loop diagrams. 

We have considered some explicit expressions following from the configuration 
space representation of water melon diagrams. In all these examples the integrals 
are ultraviolet finite. We emphasize that ultraviolet divergences add nothing new to 
the analysis. For this particular topology the structure of the ultraviolet divergences 
is particularly simple: all divergences result from the region of integration around 
small x and there are no overlapping divergences. Therefore the divergent parts of 
the diagrams can be obtained by expanding the integrand at small x and subtracting 
the corresponding singularities which simplifies the integrands drastically (making 
them effectively massless) and the integration of the divergent parts can be done 
analytically [15| within some (usually dimensional) regularization scheme. More 
complicated is to find the finite parts of the massive diagrams. Keeping this in mind 
we do not discuss the problem of ultraviolet divergences in the rest of the paper and 
concentrate on the finite parts. 

4 Reduction 

The preceding section has shown that three-loop water melon diagrams can be ex- 
pressed as configuration space integrals of a product of at most four McDonald 

15 



functions K v (x) where v need not be an integer. In this section we want to exhibit 

the three steps of how to reduce the set of necessary integrals to a smaller set. First 

we can get rid of the non-integer dimensionality of the functions by using Eq. fll7|) to 

expand into powers of e resulting in integrals containing a product of four McDonald 

functions with or without a factor ln(x). As a second step we use the relation 

11 — \ 

K n {x) = 2 tf„_i(z) + K n _ 2 {x) (39) 

x 

for McDonald functions of different orders to further reduce the integrals to integrals 
only containing KqKi, K^Kf, KoKf and Kf together with some positive powers 
of x and again with or without a factor ln(x). The last step consists in using 

^K v {x) = ~{K v ^{x) + K u+1 (x)) (40) 

for integer v (which is valid for any complex v as well) and a partial integration in 
order to reduce the necessary integrals to integrals containing only the McDonald 
function Kq(x). This reduction procedure will be considered for two different cases 
in the following two subsections. 

4.1 Reduction to integrals containing Kq(x) 

After using the recurrence relations for Bessel functions the general water melon 
diagram reduces to a linear combination of basic integrals given in Eq. (p0|). The 
subsequent reduction to even simpler integrals containing only Kq{x) and powers of 
x is done by partial integration. For this purpose we make use of the relations 

—K (x) = —K\(x) and 
dx 

^M*) = ~(K Q (x) + K 2 (x)) = -tfo(x) - ~#i(aO. (41) 

After some simple algebra for the most tedious case of four functions Ki(x) in the 
integrand we obtain the recursion relation 

d 



poo roc ri 

/ Kf(x)x q dx = - / — {Koix^Ktixfx^x 
Jo Jo dx 

r o i oo 

= - K G (x)Kl(x)x q 







/•oo />oo 

-3/ Kl(x)Kl(x)x q dx + (q-Z) K (x)KUx)x q - 1 dx. 
Jo Jo 
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The other cases are simpler and will not be written down here. For q > m, the surface 
terms of the form [KQ~ m \x) K™ (x)x q ]^' vanish. Therefore the only elements of this 
recursion are the integrals I m (q), and the recursion is expressed as 

h{q) = (q-3)I 3 (q-l)-3I 2 (q), 
h(q) = i((g-2)J 2 ( ff -l)-2/ 1 (g)), 
h{q) = l -{{q-l)h{q-l)-h{q)), 

h(q) = \qIo(q-l), (42) 
which reduces the starting integrals to our basis integrals 

roo 

I (q) = / KUx)x q dx. (43) 
Jo 

4.2 Reduction to integrals containing Kq(x)\ux 



We again use the relations (|41 ) to reduce the general water melon integrals to 
those containing only the McDonald function Kq(x). This is done again by partial 
integration. As in the previous subsection we find 

ro ° „a, , „, , f°° „„, ..dK n (x 



roo roo liKniT) 

/ KHx)x q \nxdx = - / Kf(x) — °-^x q \nxdx 
Jo Jo dx 



oo 

2/„\ t/2, 



K (x)Kt(x)x q \nx\ o -3j K£(x)Kf(x)x q \nx dx 

roo roo 

+(?-3) / K {x)K*(x)x q - 1 \nxdx+ / K (x)Kf(x)x q ~ 1 dx 
Jo Jo 

for the case with four functions K\{x) in the integrand. For integer q > m, the 
surface terms [K l ~ m ( y x)K™(x)x q lnx]g° again vanish. Therefore the recursion is 
expressed in terms of the integrals I l m {q) and is given by 

I l M = (q-3)I l 3 (q-l) + I 3 (q-l)-3I l 2 {q), 
I'M = \{{q-2)I l 2 {q-l) + I 2 {q-l)-2l[{q)) 1 

4(9) = i((ff-l)il(ff-i) + Ii(g-i)-iS(g)), 

A{q) = \{qI l Q {q-l)+h{q-l)). (44) 
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Together with Eqs. ([42]) these relations give the complete set of one-parameter re- 
currence equations for reducing a general water melon integral to a set of master 
integrals of the form Io(q) from Eq. fl43|) and 



poo 

I l ( q ) = / KUx)x q \nxdx. (45) 
Jo 



Concluding these two subsections we affirm that Eqs. fl4"5|) and (|45| ) constitute our 
basis for the evaluation of water melon diagrams with high powers of denominators. 

4.3 Efficiency of reduction 

Note that the direct reduction of a water melon diagram to the master integrals 
is rather slow within the straightforward application of recurrence relations based 
on momentum space representation. In practice the computation proceeds through 
the use of a table of integrals with given powers of the denominators. One would 
have a three-dimensional table for a given total power N if no modifications of the 
basic technique as developed in ref. || have been introduced. The number of entries 
(even when accounting for the appropriate symmetries) then grows as fast as iV 3 
which is large for the large values of N needed in some present applications. Within 
our methods one first re-expresses these integrals through a one-parameter set of 
integrals which are solved explicitly. For large N the number of entries increases as 
a first power of N (the number of elements for the Io(q) basis is given by 2[N/2] — 5 
where [z] is an integer part of z) which considerably reduces the time consumption in 
a computer evaluation. Note that in refs. [EI, different recursion techniques have 



been described which also avoid the use of the three-dimensional tables in reducing 
the water-melon diagram. For instance, in the package MATAD |21| the three-loop 
water melon diagrams are reduced to a one-dimensional table of integrals using a 
dedicated set of (momentum-based) recurrence relations. 
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5 Computation 



In this section we give some explicit formulae for computing water melon diagrams 
within the configuration space technique. The methods is best used for direct nu- 
merical computation or for constructing efficient approximation formulae for water 
melon diagrams. Therefore, as for exact analytical expressions, we mainly consider 
the D = 4 case where the initial values of the recurrence procedure can be eas- 
ily found and can be compared with numerical results. The construction of the 
£-expansion is not drastically simplified as compared to the standard momentum 
space approach and we therefore do not dwell on it here. 



5.1 N = 4 water melon 

We describe a simple way of computing an initial value of a general D = 4 (or, more 
precisely, integer dimension) water melon diagram. The quantity we need has the 
form 

I(q) = / KUx)x 2q+1 dx = I (2q + 1) (46) 
Jo 

and represents a member of our basis. 

We concentrate first on the case q = 0. Note that Eq. (ffH) is simply the result for 
a water melon diagram with a massive propagator Ko(mx) within a two-dimensional 
theory [p0|1 . The corresponding two-line water melon (master one- loop diagram) in 
momentum space has the explicit form 

n 2 (p) = 1 ln ( Vp 2 + AmI + Vp 2 ^ / 47 n 

2 r K\pp 1 \J , p 1 + 4m 2 \ v^ 2 + 4m 2 — \fp^ J 



Eq. ([46]) now becomes 

f(o)=w/a,(pw=ljfg = > < 48 » 

where we have changed the integration variable p to a new variable £ defined by the 
relation p = 2msinh(£/2) and used the standard integral p3| 

roo£«-ld£ 2 Q -1. 



o sinh £ 2 
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—T(a)((a) (49) 



with r(a) being Euler's r-function and C( a ) being Riemann's (^-function. Note that 
the two-dimensional one-loop correlator in momentum space H2 shown in Eq. 
has the simple form 

1 £ 

n 2 (2msinh(£/2)) 



Anm 2 sinh £ 

when expressed in terms of the new variable £. The integration measure becomes 
d 2 p = ATTm 2 sinh £ d£. 

Results for other values of q can be obtained by differentiating one of the two 
II2 in the integrand in Eq. ( |48|) , 

I(q) = 2nm 2 J U 2 (p)(-m 2 O p ) q Il 2 (p)d 2 p (50) 

where D p is a two-dimensional d'Alembert operator in (Euclidean) momentum 
space, O p = d 2 /dp fl dp /Ji . There is a possibility to differentiate a separate line of this 
three-loop water melon diagram that leads to different representations for higher 
moments but we find Eq. (|50| ) to be the most convenient one. We do not have a 
general analytical solution to Eq. (ISO) for arbitrary large q at the moment, although 



the solutions for some first values of q are easily available. 

Also we comment on the general structure of the basis given by I(q) in Eq. 
This basis has the form 

A q ((3) - B q (51) 

where the transcendental number £(3) is manifestly written down while A q , B q 
are rational positive numbers for any q. Several pairs of coefficients (A q , B q ) for 
q = 1,2, 3, 4, 5 are 

/ 7 3 \ / 49 27\ /63 37\ Z42777 25555\ /3101175 9304913\ 
V32'167' Vl28'64/' V32'Io7' V 2048 ' 1024 ) ' V 8192 ' 20480 / ' 

The quantity I(q) is positive for any q. Its numerical magnitude can be easily 
inferred from the asymptotic expansion of the integral given in Eq. fl46|) at large q, 

'<*> - 1 1 - ^175 +o H- (52) 

However, its analytical representation given in Eq. ([51]) reveals a rather awkward 
behaviour when analyzed numerically. Two terms with manifestly different tran- 
scendental structure taken as analytical expressions cancel each other to a large 
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extent upon numerical evaluation. Indeed, using the numerical value of the tran- 
scendental number £(3) = 1.20206 . . . , the numerical values of the pairs (A q ((3), B q ) 
for q = 1,2, 3, 4, 5 are given by 

(0.26295,0.1875), (0.460162,0.421875), (2.36655,2.3125), 
(25.1076, 24.9561), (455.052, 454.341). 

For q = 5 the numerical values of quantities of different transcendent ality coincide 
with each other up to three significant figures resulting in a cancellation on subtrac- 
tion and a large loss of precision. For larger q the cancellation of the first significant 
figures of the two numbers is more dramatic resulting in even larger losses of numer- 
ical precision. It casts some doubt on whether the presentation of results in terms 
of their manifest transcendental structure is the most preferable one. On the other 
hand the asymptotic formula Eq. ([52]) is simple and quite precise. Namely, taking 
Eq. (|5*2] ) and introducing a new parameter k which accounts for higher order terms 
in the asymptotic expansion, one has 

'<*-^-£)- 

For k = 0.97 this simple formula gives good numerical results for the basis integrals 
with an accuracy better than 1% for all q > 1. For q > 3 the relative accuracy is 
better than 10 -3 (one per mille). 

Another representation for the basis set of integrals can be obtained using the 
dispersion relation in configuration space. One obtains the relation 

K 4 (x) = (2vr) 3 / Pi (s)K (^x)ds (54) 

with p^(s) being a spectral density for the three-loop water melon diagram. Note 
that such a representation is a basis for the sum rule calculations in configuration 
space both in two-dimensional [20| and four- dimensional space-time |p4f. One ob- 



tains the following relation for the basis integrals in terms of the moments of the 
spectral density P4(s), 

/... 



I(q) = (2vr) 3 4«(g!) 2 / ^ds. (55) 
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An efficient way of computing the spectral density of water melon diagrams has 
recently been developed |25[] . However, the direct configuration space representation 
is the most convenient for a numerical evaluation. To emphasize this last remark, 
we give an asymptotic formula for the e-expansion of a water melon diagram where 
ln(x) appears in the integrand. The asymptotic formula for the log-type integrals 
reads (where we changed our notation by introducing ln(4x) for convenience) 



L\{q) = / K^{x)x 2q+l \n{Ax)dx 
Jo 



0{l/q 2 ) (56) 



9-1/2 

where ty{x) = T'(x)/T(x) is the logarithmic derivative of the T-function. 

Modifying the last term by introducing a parameter Ki as before in Eq. fl53|) we 
find that the relation 



(57) 



with Ki = 1.17 gives good results for the basis log-type integrals with an accuracy 
better than 1% for all q > 3. For q > 5 the relative accuracy is better than 10~ 3 . The 
leading order asymptotic formula presented in Eq. ( |5"7| ) is less precise for small q than 
its analog in Eq. (|53| ) in the previous case because the integrand in Eq. (p6|) is not 
positive for log-type integrals. It is also clear that an exact solution of the recurrence 
relations will be much more complicated than these simple asymptotic formulae. For 
large q the numerical cancellation of terms with different transcendental structure 
up to very high significant figures will be quite dramatic. 

The solution of the recurrence relations depends on the space-time dimensional- 
ity. In the case D = 3 the structure of the recurrence relations drastically simplifies. 
The basis set of integrals in three-dimensional space-time can be written down ex- 
plicitly. The propagator is 

D(x,m) = — e ~ mx (58) 
Atcx 

which leads to an integral basis of the form 

H (q) = (2vr) 4(A+1) / D{x,mfx q+2 dx -> — / e~ 4x x q - 2 dx = ^-T{q - 1). (59) 
Jo 4 Jo 4^ 
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The special case D = 3 can serve as a check of any general solution of the recurrence 
relations. 



5.2 N = 3 water melon: a standard sunset 

The case of a two-loop water melon (genuine sunset) is simple indeed and can be 
easily analyzed along the same lines. The corresponding basis set of configuration 
space integrals is quite analogous to the previous case and is simpler because it now 
includes only three McDonald functions, 

J n {q) = / K*- n (x)K?(x)x q dx, 
Jo 

J l (q) = / K*- n (x)K?(x)x q \nxdx. (60) 
Jo 

The reduction to the basis set of integrals analogous to the case of three-loop water 
melons given in Eqs. fl4"3"| ) and ([P|) can be now readily obtained. 

The basic initial integral (the basic sunset Bg) for the recurrence relation has 
the explicit form 

with n~2(p) given by Eq. (f47j). Here M is a mass of the third line which is kept 
different from the other two with masses m. By differentiating with respect to 
M, any positive power of propagator (and/or power of x in configuration space 
representation) can be obtained. After changing the integration variable as in the 
preceding subsection one finds an explicit representation 

roa 

Bs = Jo Am? sinh 2 (£/2) + M 2 ' ^ 
After a change of variables the integration can be done and can be reduced to a 
polylogarithm function. Namely, for t = one has 

B r zjz _ 1 r 1 ^ tdt m 

S Jo 4m 2 sinh 2 (£/2) + M 2 m 2 Jo 1-2-ft + t 2 1 ' 

where 7 = 1 — M 2 /2m 2 . The last integral is evaluated to give 

Intdt Li 2 (l/ti) - Li 2 (l/t 2 ) 



o 1 - 27* + t 2 tx- t 2 
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(64) 



where t lf2 = 7 ± a/7 2 — 1 and Li 2 (z) is the dilogarithm function, 

00 k 

^2(2) = £73, 1*1 <!■ 
k=i K 

The differentiation with respect to M is now straightforward and can be performed 
with a symbolic manipulation program. 

In the case M = 2m the integration simplifies because the two independent 
parameters M and 2m on which the integrand depends, coincide (degenerate case). 
The integral is then reduced to a special case of Eq. fl4"9|) . We have 7 = — 1 and 

- 1 \ntdt 



(1+t) 2 



In 2 (65) 



which leads to 

B S (M = 2m) = (66) 
m A 

For the case M = m the standard result - Clausen's polylogarithm Cl2(vr/3) - is 
reproduced (see e.g. ref. ||). Indeed, 7 = 1/2 and — exp(±i7r/3). Eq. ([64]) now 
becomes 

r 1 Intdt 2 T T , . /q . , . 

and with using the definition of Clausen's polylogarithm, C^(^) = ImLi2(e ) one 
finds 

Bs {M = m) = -^=C1 2 (I) . (68) 



m 



v/3 V3 



5.3 Generalization to the spectacle topology 

In this subsection we give a formula for a more general topology when only one 
propagator is removed from the initial three-loop bubble diagram. In the original 
classification of ref. these are class E diagrams belonging to the spectacle topol- 
ogy. The formula obtained in this subsection is efficient for numerical integration 
though we did not find any analytical solution. The main obstacle of generalizing 
the configuration space technique to a general multi-loop diagram is the angular 
integration. The configuration space technique proved to be rather successful for 
general diagrams in the massless case but it brings no essential simplification in 
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the general massive case (see e.g. ref. [[27j). However, for special configurations the 
angular integration can be explicitly performed with a reasonably simple integrand 
left for the radial integration. The diagrams of spectacle topology give an example 
of such a configuration. 

The configuration space expression of a spectacle topology diagram written in a 
form suitable for our purpose is (see Fig. 2(a)) 

/ D(x — y,m)D(x,m ) 2 d D xD(y,m) 2 d D y. (69) 

The key relation for a drastic simplification of the configuration space integral with 
the spectacle topology is the addition theorem for Bessel functions allowing one 
to perform some angular integration explicitly. One needs to integrate over the 
relative angle in the propagator D(x — y,m). In the handbook of Gradshteyn and 
Rhyshik |T7f one finds 



« = 2-m-TM J> + k) ^^ Z " +k[ r ] CU^V) (70) 
n k=0 P r 

where C% are the Gegenbauer polynomials, Z is any of the Bessel functions J, Y, 

H W or #(2) 



R = \jr 2 + p 2 — 2rp cos tp (71) 

and r > p. For r < p the arguments of the Bessel functions on the right hand side 
of Eq. ([70) should be interchanged. Writing R — \r — p\ we have for r > p 

Z v {m\r - p\) " J u+k {mp) Z u+k {mr) 

— : — = 2m T{u)2^{u + k) C fc (cos^). 72 

\r - p\ u ^ p v r v 

Using this relation for Z = and substituting m = e tn ^ 2 for the purpose of 

analytic continuation in order to obtain a relation for the modified Bessel functions 

K and /, we find 



K ^ r ~J ] = 2^F(A) t(X + k) h+ ^ p) Kx : k x (r) Cjiccs (73) 
where we have 



| r _p|A p\ r x 



71 



Kx{z) = je^H?(iz), 
I x ( z ) = e - nXl/2 J x {e wi/2 z) for -7i < aigz < 
I x ( z ) = e 3nXl/2 J x (e' 3m / 2 z) for ^<aigz<n. (74) 
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Using the orthogonality relations for Gegenbauer polynomials (see the Appen- 
dix), the sum disappears after integration over the relative angle and only one term 
contributes. We obtain 



- 2T (^—x —f(ITTj C " {1) 

= (^f^ r> P , (75) 

where the first equality is a consequence of the orthogonality relation with the trivial 
factor Cq(1) = 1. This result allows one to write down an expression for any 
spectacle-type diagram in the form of a two-fold integral with a simple integration 
measure 



poo roo 

/ D(x,m) 2 x 2X+1 dx / D(y : m) 2 y 2X+1 dy 
Jo Jo 



X 



where 0(x) is the standard step-function distribution. 

Note that the integration measure D(x,m) 2 x 2X+1 dx allows one to perform the 
integration by using efficient integration routines for a numerical evaluation. The 
form of the weight function is close to e~ ax x a which makes the use of Laguerre 
polynomials a convenient choice within the Gaussian numerical integration method. 
Any modified propagator (with any power of the denominator) can be used as a 
factor in the integration measure D(x,m) 2 x 2X+1 dx which makes this representation 
universal and useful for the case of high powers of denominators of the lines associ- 
ated with pairs (x, 0) and (y, 0) of space-time points. If the angular structure of the 
diagram is preserved, the generalization to higher loops in the expressions for the 
radial measures is straightforward. 

As an illustration of this technique we present an example of the evaluation of 
a spectacle diagram. Consider an integer dimension space-time which, without loss 
of generality, can be chosen to be two-dimensional (an odd number of dimensions is 
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trivial because the propagators degenerate to simple exponentials). The spectacle- 
type three-loop diagram can be obtained in a closed form. Indeed, in the momentum 
space representation we have 

s m=jM& <77) 

for the basic spectacle diagram S with n 2 Qo) taken from Eq. ([|7|) and the mass M of 
the connecting propagator kept different. After the substitutions p = 2msinh(£/2) 
and t = we have 

^ (M) = 2^L {l-t*)£-tyt + &)' (?8) 
Performing the integration we finally obtain 

S(M) = /(tl) ~ /(t2) (79) 

t\ — t-2 

with 

= 8tLi 3 (l/t)-(t + 7)C(3) 
JU 87rm 4 (t 2 -l) 

where t\p = 7 ± y/'y' 2 — 1 with 7 = 1 — M 2 /2m 2 as before. Li 3 (z) is the trilogarithm 
function 

oo 

Li 3W = Ep. 1^1 < L 

k=l K 

For M = 2m the integral in Eq. (fTTf ) simplifies as in the case of the sunset diagram 
and one finds a simple answer in terms of the standard (in the present context) 
transcendental numbers In 2 and C(3), 

For the actual value of the mass M = m we obtain a result including the next 
Clausen polylogarithm Ci3(27r/3). As one can conclude from this expression, the 
conjugate pair of the sixth order roots of unity, exp(±27ri/3) plays an important 
role in this case again in accordance with the general analysis of ref . || . The origin 
of the appearance of the the sixth order roots of unity as the parameters of the 
analytical expressions of the diagrams lies in the mismatch of masses along the lines 
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of the diagrams. However, the exceptional case M = 2m, where one line has the 
double mass of the other lines (which results in the drastic simplification) also keeps 
us within the set of the sixth order roots of unity. The key parameter in this case is 
simply the natural number 1 which definitely is one of the sixth order roots of unity. 
Turning to the configuration space representation we find 



S(M) = / xK (mx) dx / yK (my) dy x 
Jo Jo 

\K Q (Mx)h{My)e{x -y) + K (My)I (Mx)6(y - x)) (80) 



for the basic spectacle diagram. An explicit numerical integration of Eq. (|80j) shows 
coincidence with the analytical result in Eq. ([F^) which we checked for arbitrary 
values of M. In this example the analytical result has a rather simple form which 
is not true if high powers of denominators enter. Then the corresponding one-loop 
insertions are rather cumbersome and an explicit integration in configuration space 
is more convenient. 

As a last topic of this subsection we demonstrate how the square of Clausen's 
polylogarithm Cl2(vr/3) 2 can emerge at the level of spectacle topology diagrams. The 
transcendental number Cl2(vr/3) 2 characterizes the analytical results for three-loop 
bubbles. Its presence was discovered in the impressive treatise of David Broadhurst 
on the role of the sixth order roots of unity for the transcendentality structure of 
results for Feynman diagrams in quantum field theory ||. 

We consider the spectacle diagram in the form shown in Fig. 2(b). Here the 
expression for the generalized middle line is a product of the one-loop propagator of 
Eq. ( fi?D and the standard particle propagator. We express the one-loop propagator 
in Eq. p7| ) using a dispersion representation with the spectral density P2{s) and 
obtain 

. , _ , , ~ , . 1 1 f°° p 2 (s)ds 

Loop(p,m) x Propagator (p, m) = n 2 (p)^— - = ^— / ■ — - 

p 2 + m 2 p 2 + m 2 Jim? s + p 2 

= r^f ' ' ). (si, 

Jim 2 s — m z \p 2 + m z s + p z J 
Taking only the first term (which is sufficient for obtaining the result we are looking 
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for) one has 

/= i reM* (82) 

pZ _|_ m Z J 4m 2 g _ m Z 

which leads to the sunset diagram after the two other line shown in Fig. 2(b) have 
been added with a normalization factor given by the integral. One factor Cl2(7r/3) 
results from integrating the overall sunset diagram which is composed of the prop- 
agator (p 2 + m 2 ) -1 from Eq. fl32|) with the two other lines of the diagram shown in 
Fig. 2(b). The second factor Ci2(7r/3) has to be found in the normalization factor 
given by the integral in Eq. (0). Note that the very structure of this contribution 
- the square of a number which first appeared at the lower loop level - suggests a 
hint for its search. It should emerge as an iteration of a lower order contribution in 
accordance with the iterative structure of the i?-operation (see e.g. ref. [^8|) which 
forms a general framework for the analysis of multiloop diagrams. The following 
consideration confirms this guess. Consider the quantity 

N= f°° P2(s)ds 
Jim 2 s — m? 

and take p2(s) to be the spectral density in D-dimensional space-time (see e.g. 
Eq. (65) of ref. 

P2(S) = 2^+1/^ + 1/2)^' ^ > ^ (83) 

Now consider a first order contribution of the expansion in e near the space-time 
dimension D = 2. The expansion in A = — e results in 

(s - 4m 2 )- E ~ 1 /2 1 / , (s-Am 2 \ 2 A . . 
ITr = / 1-gln — +Oe 2 , 84 

so the relevant first order term in e is 

. . ln((s - 4m 2 )/m 2 ) . K . 

A £ p 2 {s) = vv ;/ ; (85) 

2-7T Js(s - Am 2 ) 



where p = m has been chosen for convenience. Now we change the variable according 
to 

= 2mcosh(£/2), t = e" f (86) 
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to obtain 



, . ., (Int-2In(l -t))t 

A«„ 2 (W co 8 h 2 ( f /2)) = 2mni(1 [ _ t2) " ■ (87) 



For the quantity in question we find 

A N = [°° A ^ ds - _L_ f 1 (lnt-21n(l-t))dt 

Am* s-m 2 ~~ lum 2 Jo 1+t + t 2 ' 1 ' 

The roots of the denominator of the integrand in Eq. fl88l) are now £3^ = exp(±27rz/3) 
which again is a conjugate pair of the sixth order roots of unity. After integrating 
this equation we readily find 

A ' N = dvl( Im ( Li2 ( c2 "* /3 )-^(^" /3 ))-i'«3) 

bHIK"4 w 



7rm 



Using the relation 



we finally obtain 



*'f)-b-(f 



A ^-^b( ci '©-i in3 )- (9o » 

Therefore, in the first order of the e expansion of the spectacle diagram we found 
this remarkable contribution proportional to Cl2(vr/3) 2 . In our calculation it emerges 
naturally as the iteration of the lower order term. Originally this contribution had 
been guessed and confirmed in ref. || by a direct numerical computation of the 
finite part of the general three-loop bubble in four- dimensional space-time. 

Note also that besides and Cl2(7r/3) 2 there is one more nontrivial transcenden- 
tal number necessary for the analytical computation of the finite part of the general 
three-loop bubble in four-dimensional space-time as has been proven in ref. || . This 
number was chosen in ref. as a two dimensional nested sum of the form 

^ (-l) m cos(27m/3) 

^3,1 - 2^ 3 



m >n>0 17171 



The analogous contribution may also appear in the first order e-expansion of the 
spectacle diagram. We hope to devote a separate publication to this point. 
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Therefore there are solid arguments that all basic transcendental numbers nec- 
essary for the analytical computation of three-loop massive bubbles can be found 
at the level of much simpler topology than the general tetrahedron. One may need 
only spectacle diagrams as the largest set to find all necessary transcendental num- 
bers. These discoveries lead us to a conjecture about the extension of recurrence 
relations for three-loop bubbles beyond those found within the integration by part 
technique. Indeed, if all results can be obtained within the spectacle topology only, 
it looks plausible that there is a procedure that can perform a reduction of the gen- 
eral tetrahedron diagram to a simpler set where the spectacle topology is the most 
complicated one. The manifest form of this procedure is not known at present. 

6 Conclusion 

We have formulated a new representation for some of the massive three-loop vacuum 
diagrams with the simple water melon topology. We have obtained numerical values 
for some master integrals within our representation. The computation of diagrams 
with large powers of propagators is reduced to a linear (one-parameter) set of basis 
integrals. Our representation is simple and provides a tool for an efficient numerical 
evaluation of such diagrams. We have given asymptotic estimates which provide 
an accuracy better than one per mille for all basis integrals with q > 3 in the 
case of the leading order of e-expansion and with q > 5 in the case of the first 
order of e-expansion. The generalization to higher loops is straightforward. We 
gave strong arguments that all nontrivial transcendental numbers necessary for the 
computation of three-loop bubbles (which were identified in ref. || by using a high 
precision numerical integration) already appears at the level of spectacle topology 
diagrams. 

As a by-product of our analysis we obtained analytical results for a number of 
definite integrals involving products of four McDonald functions of different orders 
with powers and logarithms (Sec. 3). 
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Appendix 

We give some formulae for dealing with Gegenbauer polynomials in Sec. 5.4 that 
can be found e.g. in refs. 2S|. The Gegenbauer polynomials obey the orthog- 



onality relations 



J C*(x! ■ x 2 )C^(x 2 ■ x 3 )dx 2 = ^ 5 mn C*{xi -x 3 ) (J dx 2 = 1) 
or, written in another form, 

f 27T A+1 A 

jciia- b)C} 2 (b ■ c)dn b = 8 jlh nx + 1)ji + ) C*(a ■ c). 
In particular, we have Cq(x) = 1, C*(x) = 2Xx and 

(j + l)C} +1 {x) = 2(j + X)xC}(x) - (j + 2A - VC^x) 

together with 



tA m _ r(y + 2A) 
3 1 ] j!F(2A) 



The characteristic polynomial is given by 

oo 

(t 2 -2tx + \y x = Y J t 3 C^(x). 



3=0 
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